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The  mechanical  properties  of  a  cross-linked  polydimethylsiloxane  (PDMS)  network  swollen  with  nonreactive  entangled 
PDMS  solvent  was  previously  studied  experimentally.  In  this  article,  we  use  the  discrete  slip-link  model  to  predict  its 
linear  and  nonlinear  rheology.  Model  parameters  are  obtained  from  the  dynamic  modulus  data  of  pure  solvent.  Network 
rheology  predictions  also  require  an  estimate  of  the  fraction  and  architecture  of  dangling  or  inactive  strands  in  the  net¬ 
work,  which  is  not  directly  measurable.  The  active  strand  fraction  is  estimated  from  dynamic  modulus  measurements, 
and  the  molecular  weight  is  adjusted  to  fit  the  dynamic  modulus  data.  Then,  the  nonlinear  rheology  can  be  predicted 
without  adjustments.  These  successful  predictions  strongly  suggest  that  the  observed  rheological  modification  in  the 
swollen  blend  arises  from  the  constraint  dynamics  between  the  network  chains  and  the  dangling  ends.  ©  2014  American 
Institute  of  Chemical  Engineers  AlChE  ./,  60:  1372-1380,  2014 

We  dedicate  this  article  to  Prof.  R.  Byron  Bird  on  the  occasion  of  his  90th  birthday.  We  attempt  to  follow  the  twin 
exhortations  of  Bob  to  be  both  mathematically  rigorous  and  practically  relevant. 
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Introduction 

Cross-linked  polymer  networks  swollen  with  polymeric 
solvent  have  shown  adaptive  mechanical  response.  This  fre¬ 
quency  dependent  response  makes  these  gels  desirable  for 
applications  such  as  tissue  simulants  for  ballistic  testing.1 
Recently,  the  dependence  of  solvent  molecular  weight  and 
concentration  on  the  mechanical  and  rheological  properties 
of  an  end-linked  polydimethylsiloxane  (PDMS)  gel  that  con¬ 
tains  PDMS  soluble  material  were  studied  by  Mrozek  et  al.2 
In  that  work,  the  focus  was  on  an  end-functionalized  PDMS 
precursor  that  is  much  larger  than  the  entanglement  molecu- 
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lar  weight,  resulting  in  a  low  chemical-cross-link  density  and 
a  modulus  that  is  determined  by  the  degree  of  entanglement. 
As  a  result,  a  change  in  the  entanglement  screening  or 
enhanced  looped  formation  will  have  a  significant  impact  on 
the  mechanical  properties  of  the  gel.  The  networks  are 
formed  in  the  presence  of  nonreactive  PDMS  solvents  with 
molecular  weights  broadly  spanning  the  molecular  weight  of 
entanglement  for  PDMS.  The  increased  frequency  depend¬ 
ence  of  the  mechanical  response  and  improved  gel  fracture 
toughness  with  entangled  solvent  is  expected  to  benefit  many 
applications  including  electronic  device  encapsulation,  pres¬ 
sure  sensitive  adhesives,  and  tissue  surrogates  for  biomedical 
engineering. 

Gels  swollen  with  high-molecular-weight  solvent  were 
previously  modeled  using  standard  polymer  coarse-grained, 
bead-spring  “Kremer-Grest”  models.3-5  These  simulations 
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Figure  1.  Sketch  of  how  the  swollen  PDMS  network  is  modeled  using  the  DSM. 

In  the  DSM,  each  component  is  simulated  independently  using  a  self-consistent  mean  field  that  accounts  for  the  effect  of  other 
types  of  chains.  The  components  of  the  swollen  network  modeled  here  are,  A:  SC,  B:  DC,  C:  NC.  The  circles  represent  binary 
entanglements  and  the  crosses  represent  cross-links.  Both  of  which  are  fixed  in  space  for  Green-Kubo  calculations  or  moved 
affinely  for  flow. 


have  proven  to  be  a  useful  technique  to  study  microscopic 
topology  and  dynamics  of  lightly  to  moderately  entangled 
polymer  systems.  However,  the  models  have  limitations 
when  studying  very  high-molecular-weight  polymer  chains 
where  entanglement  dynamics  become  dominant,  because 
the  time  scales  of  interest  are  frequently  longer  than  what  is 
reachable  computationally.  Attempts  have  also  been  made  to 
apply  tube  models6'7  to  describe  the  elasticity  of  both  dry8’9 
and  swollen10'11  polymer  networks.  However,  the  time- 
dependent  entanglement  effect  on  the  relaxation  of  poly¬ 
meric  networks  remains  unresolved  in  those  theories.12 

In  this  article,  we  apply  the  discrete  slip-link  model 
(DSM)  introduced  by  Schieber  and  coworkers, 13-17  which 
was  developed  to  describe  the  dynamics  of  flexible  entangled 
polymer  liquids  and  gels.  The  model  is  able  to  predict  the 
linear  viscoelasticity  of  monodisperse  linear,  polydisperse 
linear,  and  branched  systems.  Recently,  the  model  also 
showed  mostly  good  agreement  with  dielectric  relaxation 
experiments.18  Without  any  adjustment,  the  DSM  shows 
excellent  agreement  with  shear  flow  experiments  and  elonga- 
tional  flows  with  elongation  ratios  of  up  to  10-20,  as  shown 
recently.19 

The  model  was  first  applied  to  cross-linked  polymer  gels 
with  entanglement  dominated  dynamics  by  Jensen  and 
coworkers.211  In  that  work,  the  effect  of  network  imperfections 
(dangling  network  chains)  on  the  dynamic  modulus  of  the  gels 
was  systematically  studied  using  the  DSM. 

In  this  article,  we  expand  on  that  work  and  use  the  DSM 
to  estimate  the  composition  and  predict  the  relaxation  and 
nonlinear  mechanical  response  of  cross-linked  polymer  net¬ 
works  swollen  with  polymeric  solvent.  The  main  objective  is 
to  test  a  methodology  for  predicting  and  designing  the  rheo¬ 
logical  properties  of  gels.  Toward  this  end,  we  test  the  feasi¬ 
bility  of  using  the  model  to  characterize  the  fraction  of  gel 
strands  that  are  active  from  the  dynamic  modulus,  and  check 
if  the  nonlinear  rheology  can  then  be  predicted.  The  experi¬ 
mental  system  modeled  is  a  cross-linked  PDMS  network  in 
the  presence  of  high-molecular  weight  PDMS  linear  solvent. 
When  accounting  for  imperfections  in  the  cross-linked  net¬ 
work,  this  application  requires  the  coupling  of  three  catego¬ 
ries  of  architectures:  (1)  entangled  network  chains  (NC);  (2) 
dangling  network  chains  (DC);  and  (3)  entangled  solvent 
chains  (SC).  Additionally,  we  account  roughly  for  the  poly- 
dispersity  in  the  molecular  weight  distribution  of  each  of  the 
components  in  the  network  blend. 

The  article  is  organized  as  follows,  in  Model:  DSM  for 
Blends  Section,  we  summarize  the  main  aspects  of  the  DSM 


with  an  emphasis  on  how  the  interaction  between  different 
types  of  chain  architectures  and  molecular  weights  are 
handled  in  the  single-chain  self-consistent  mean-field 
approach.  In  Blend  Composition  Section,  the  details  about 
the  composition  of  the  system  being  modeled  are  presented. 
There,  we  describe  the  way  in  which  we  account  for  the  pol- 
ydispersity  of  each  component  in  the  network  blend.  Also, 
the  methodology  to  estimate  the  fraction  of  imperfections  in 
the  network  structure  (dangling  chains)  is  described.  In 
Determination  of  the  DSM  Parameters  for  PDMS  Section, 
the  two  parameters  that  are  required  in  the  DSM  to  predict 
rheological  properties  are  estimated  using  dynamic  modulus 
data  for  polydisperse  entangled  PDMS  liquid.  These  same 
parameters  can  then  be  used  to  model  all  other  PDMS  chain 
architectures  and  molecular  weights.  In  Dynamic  Modulus 
Predictions  of  the  Swollen  Network  Section,  the  dynamic 
modulus  of  a  cross-linked  PDMS  network  swollen  with  poly¬ 
meric  PDMS  solvent  is  predicted  using  the  DSM  and  com¬ 
pared  to  experimental  data  by  adjusting  the  only  remaining 
parameter:  the  molecular  weight  of  the  dangling  network 
chains.  An  analysis  on  the  effect  of  high-molecular-weight 
solvent  in  the  linear  viscoelastic  response  of  cross-linked 
polymer  networks  is  given.  In  Predictions  of  Simple  Elonga¬ 
tion  of  Swollen  Network  Section,  a  simple  elongation  experi¬ 
ment  of  the  swollen  PDMS  network  is  tested  against  the 
DSM  prediction.  Moreover,  important  insight  into  the  spe¬ 
cific  dynamics  of  each  component  in  this  complex  blend  is 
gained. 

Model:  DSM  for  Blends 

The  DSM  is  a  single-chain  mean-field  model  for 
entanglement-dominated  polymer  dynamics  proposed  by 
Schieber  et  al.14  and  evolved  by  Khaliullin  and  Schieber.16 
After  a  brief  outline,  we  apply  the  model  to  polymeric  gels 
swollen  with  high-molecular-weight  solvent. 

The  network  and  solvent  chains  are  described  by  random 
walk  statistics.  This  is  expected  to  hold  for  polymeric 
chains  with  contour  length  and  entanglement  spacing  larger 
than  several  Kuhn  steps.  The  model  is  a  single-chain  mean- 
field  model,  where  Figure  1  shows  cartoons  of  the  coarse- 
graining  process  of  the  real  network  and  solvent  chains. 
The  chains  are  first  coarse-grained  by  random  walk  statis¬ 
tics;  then,  they  are  further  coarse-grained  to  the  primitive 
path  defined  by  the  entanglements,  assuming  that  the  chain 
relaxation  is  slower  than  the  relaxation  of  a  strand  between 
entanglements.  The  entanglements  are  randomly  distributed 
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along  the  chain  with  uniform  probability  1/(1+/?),  where 
/?  is  a  model  parameter  related  to  the  entanglement  density. 
The  primitive  path  defines  the  entanglement  spacing.  The 
number  of  Kuhn  steps,  Nh  in  the  ith  strand  fluctuates  due 
to  Kuhn  step  shuffling  through  the  slip-links,  a  process 
called  “sliding  dynamics”  (SD).  When  this  process  occurs 
at  the  ends  of  a  solvent  chain,  or  at  the  uncross-linked  end 
of  a  dangling  chain,  entanglements  can  be  destroyed  or 
created. 

At  equilibrium,  the  probability  density  for  a  given  chain 
conformation  is  governed  by  the  modified  Maxwell-Boltz- 
mann  relation,  p]eq  (ST)  where  y  denotes  chain  length  and 
=  (Z,  {Nj},  {Q/}7  {rpD  })  determines  the  chain  conformation 
in  the  level  of  description  of  the  DSM.  Where  Z  is  the  num¬ 
ber  of  entangled  segments,  and  Q,  is  the  vector  connecting 
entanglements  i  —  1  and  i.  pyeq  (Q)  is  given  by 


Pjq(0)  = 


exp 


F(£lj 


kBT 


Up 


,CDfTCO)  (1) 


where  /=(  1  +  1//?)Wk  1  is  the  normalization  constant,  8(i,f) 
—Sjj  is  the  Kronecker  delta  function  which  is  responsible  for 
the  conservation  of  Kuhn  steps  in  a  chain,  N]K  is  the  total 
number  of  Kuhn  steps  in  a  chain  of  type,  y,F(Q)  is  the  free 
energy  of  a  chain  with  conformation  D,  pCD  (tcd  )  is  the  con¬ 
straint  dynamics  (CD)  probability  density  of  entanglement 
lifetimes,  which  in  the  mean-field  description  accounts  for 
the  creation  and  destruction  of  entanglements  anywhere 
along  the  chain,  kB  is  Boltzmann’s  constant,  and  T  is  the 
absolute  temperature. 

The  random  walk  of  N,  Kuhn  steps  between  two  slip-links 
of  separation  \Qj\=Qi  can  be  approximated  by  a  Gaussian 
free  energy  derived  by  Schieber  et  al.14 

Fs(Qr,N,)_  3 Q]  |  3 

kBT  2 N,a\  2 


2nNja^ 


(2) 


where  aK  is  the  length  of  a  Kuhn  step.  The  free  energy  for  a 
chain  is  the  sum  of  the  free  energies  of  the  segments  in  the 
chain 


Z— a 

F(n)=EV.(G/,ty).  (3) 

i—  1 

where  a  takes  a  value  of  0  for  a  network  chain,  1  for  a  dan¬ 
gling  chain,  and  2  for  a  solvent  chain.  For  a  network  chain, 
both  strand  ends  are  fixed  and  accounted  for  in  the  free 
energy  of  the  chain.  For  the  dangling  chain,  only  one  end 
strand  is  fixed.  And  for  the  solvent  chain,  neither  end  strand 
contributes  to  the  chain  free  energy.  This  is  because  the  end 
strand  samples  all  possible  orientations,  averaging  to  a  zero 
free  energy.  The  free  energy  of  an  entangled  strand,  Eq.  2,  is 
valid  for  polymer  chains  which  are  not  stretched  beyond 
approximately  1/3  of  their  contour  length. 

The  entanglement  activity  of  the  surrounding  chains  is 
responsible  for  fluctuations  in  the  number  of  entanglements 
on  the  chain.  At  equilibrium,  Kuhn-step  shuffling  between 
entangled  strands  is  due  to  Brownian  forces  and  free  energy 
differences. 

The  characteristic  lifetime  of  the  ith  entanglement,  tpD ,  is 
introduced  to  implement  a  mean-field  self-consistent  realiza¬ 
tion  with  independent  chains  in  the  ensemble.  The  lifetimes 
rpD  are  chosen  from  the  distribution  of  lifetimes,  pCD .  This 


distribution  is  determined  self-consistently  from  destruction 
of  entanglements  in  a  system  by  SD.  Entanglements  between 
two  entangled  network  chains  are  permanently  trapped; 
hence,  rpD  is  infinitely  large  and  for  those  Pnc(tcd)  = 
2i5(1/tcd)/(tcd)2.  The  pCD  for  the  swollen  network  blend 
considered  in  this  work  is  given  by 

Npc  Nsc 

„CD  /_CD  \ _  DC  CD  .  ,  Y^^SC  „CD 

P  (T  )-WDC  2^Wy  PDC.y+wSC  )  Psc  ,y 

7=1  y=l 

+wncEmv,CPnc,t  (4) 

y=l 


where  N  is  the  number  of  components  in  the  polydisperse 
distribution  of  molecular  weights  of  each  type  of  architec¬ 
ture  in  the  blend.  Destruction  and  creation  of  entanglements 
by  CD  occur  anywhere  along  the  chain,  which  is  different 
from  creation  and  destruction  of  entanglements  by  SD  which 
can  only  occur  in  dangling  ends.  Therefore,  in  network 
chains  in  a  blend,  such  as  the  one  considered  here,  entangle¬ 
ments  can  be  created  and  destroyed  by  CD.  In  the  dangling 
chains  of  the  network  (network  imperfections),  entangle¬ 
ments  can  be  created  and  destroyed  by  SD  only  in  the 
uncross-linked  end  and  by  CD  anywhere  along  the  chain. 
For  the  solvent  chains,  entanglements  are  created  and 
destroyed  by  SD  on  both  ends  of  the  chain  and  by  CD  any¬ 
where  along  the  chain. 

Therefore,  in  the  DSM  each  component  of  the  swollen 
network  blend  is  modeled  in  a  self-consistent  realization 
with  independent  chains  in  the  ensemble.  The  effect  of  other 
chains  is  given  by  the  self-consistent  mean-field  defined  in 
Eq.  4.  This  concept  is  illustrated  in  Figure  1.  In  other  words, 
in  the  DSM,  each  component  of  the  complex  swollen  net¬ 
work  blend,  with  different  architectures  and  polydispersity  in 
each  of  the  components,  can  be  realized  independently  while 
the  effect  of  other  architectures  and  molecular  weights  is 
given  by  the  CD  mean-field. 

The  probability  density  in  Eq.  1  evolves  in  time  according 
to  a  differential  Chapman-Kolmogorov  equation.13  Two 
types  of  calculations  can  be  performed,  equilibrium  (or 
Green-Kubo)  calculations  in  which  the  rate  of  deformation 
tensor21'22  is  set  to  zero  and  the  autocorrelation  function  of 
stress  at  equilibrium  is  followed;  or  flow  calculations  in 
which  a  specific  flow  field  is  applied  and  the  stress  as  a 
function  of  time  is  recorded. 

The  stress  tensor  is  obtained  from  thermodynamics23 


N 


where  a  has  value  0  for  a  network  chain,  1  for  a  dangling 
chain,  and  2  for  a  solvent  chain.  After  calculating  pCD(rCD) 
and  using  it  in  Eq.  1,  the  relaxation  modulus  of  the  blend, 
G(t)  is  predicted  using  the  Green-Kubo  expression 


G7M=— rET«(°)Txy(?))eq,7  (6) 

where  (...)  is  an  ensemble  average,  Txy(t)  can  be  any  off- 
diagonal  stress  tensor  component,  and  ny  is  the  number  of 
chains  per  unit  volume.  Note  that  the  stress  tensor  of  the  sys¬ 
tem  is  a  sum  of  the  stresses  of  its  individual  chains  per  unit 
volume  as  can  be  seen  in  Eq.  5.  From  this,  it  follows  that 
the  relaxation  modulus  of  the  blend  can  be  expressed  in 
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Table  1.  Molecular  Characteristics  and  Parameters  for  the 
GEX  of  the  Nonreactive  PDMS  Soluble  Solvent  Chains 
(m-PDMS)  and  Gel  Precursor  Chains  (v-PDMS) 


Name 

Mw 

(kDa) 

Mn 

(kDa) 

Mw/Mn 

(kDa) 

mp 

a 

b 

m-PDMS 

187.5 

116.8 

1.61 

263.7 

35074.1 

2.057 

0.810 

v-PDMS 

116.1 

74.8 

1.55 

165.7 

1984.09 

3.535 

0.536 

terms  of  the  relaxation  modulus  of  each  of  the  components 
in  the  blend  as 

Npc  A^sc 

G{t)  =wDC  ^2  w°c  G°c  (t)  +wsc  G*c  (f) 

y— 1  y=l 

A'nc 

+  H-NC  J>*cGf  (t)  (7) 

y=t 


Blend  Composition 

One  of  the  main  challenges  is  to  model  a  blend  with  com¬ 
plex  and  partially  unknown  composition.  More  specifically, 
here,  we  consider  a  realistic  network  with  imperfections. 
These  imperfections  are  caused  by  incomplete  cross-linking 
during  curing  or  structures  that  are  not  active.  Additionally, 
most  commercially  available  precursors  and  high-molecular- 
weight  solvents  have  significant  polydispersity  in  their 
molecular  weight  distributions.  In  this  section,  we  propose  a 
procedure  to  quantify  approximately  these  two  compositional 
parameters  theoretically. 

The  specific  system  to  be  modeled  was  studied  experimen¬ 
tally  by  Mrozek  et  al.2  and  is  composed  of  a  gel  formed  by  a 
vinyl-terminated  v-PDMS  precursor  that  is  much  larger  than 
the  entanglement  molecular  weight  resulting  in  a  low  chemical 
cross-link  density  and  a  modulus  that  is  determined  by  the 
degree  of  entanglement.  The  networks  are  formed  in  the  pres¬ 
ence  of  nonreactive  PDMS  methyl-terminated  m-PDMS  sol¬ 
vents  broadly  spanning  the  molecular  weight  of  entanglement 
for  PDMS.  The  samples  are  cured  for  72  h  at  82°C  in  the  pres¬ 
ence  of  a  platinum-cylcovinylmethylsiloxane  complex. 

Our  goal  is  to  predict  simultaneously  the  dynamic  modulus 
and  elongation  experiments  of  one  of  the  swollen  network 
blends  studied  experimentally  by  Mrozek  et  al.2  In  that  article, 
PDMS  networks  with  different  amounts  of  high-molecular- 
weight  PDMS  solvent  were  prepared.  Here,  we  focus  on  the 
blend  with  50%  loading  of  high-molecular-weight  solvent. 


Polydispersity  of  components 

The  molecular  weight  distribution  (MWD)  of  the  precur¬ 
sor  chains  and  the  PDMS  solvent  chains  can  be  measured 
experimentally  using  gel  permeation  chromatography.  For 
the  specific  system  studied  in  this  article,  the  results  of  these 
measurements  are  shown  in  Table  1.  The  molecular  weight 
distribution  obtained  from  these  experiments  can  be  fit  to  a 
generalized  exponential  function  (GEX)24 


W(M)  = 


b 


(8) 


where  the  parameters,  mp,  a,  and  b,  listed  in  Table  1,  are 

r(2±i) 

related  to  the  moments:  M^=mp  ,Mw=mp  , Mz  = 

r\b)  r(— ) 


eM 


Note  that  the  solvent  is  long  enough  to  entangle,  as  Mt 
Qi  29  kDa  for  PDMS  chains.2  The  MWD  for  the  v-PDMS 
precursors  are  shown  in  Figure  2B  and  the  MWD  for  the  m- 
PDMS  chains  are  shown  in  Figure  2A. 

It  is  convenient  to  approximate  the  GEX  distribution  with 
a  set  of  discrete  molecular  weights  for  calculating  rheology. 
The  simplest  discrete  approximation  is  bidisperse,  where  the 
first  and  second  moments  were  matched.  We  find  this  simple 
approximation  to  be  sufficient  to  reproduce  experimental 
results.  In  principle,  one  can  use  a  discrete  approximation  of 
the  GEX  function  with  as  many  molecular  weights  as  neces¬ 
sary  to  obtain  a  good  prediction  of  experimental  results. 
However,  this  also  increases  the  computational  cost.  The 
bidisperse  approximations  are  plotted  together  with  the  GEX 
fit  of  the  molecular  weight  distribution  for  the  solvent  and 
precursor  chains  in  Figures  2A,  B,  respectively. 

The  molecular  weight  distribution  of  the  dangling  network 
chains  must  also  be  estimated,  which  cannot  be  measured 
experimentally,  therefore,  we  model  them  as  monodisperse 
chains  with  a  molecular  weight  chosen  to  obtain  good 


Figure  2.  (A)  Bidisperse  approximation  of  the  molecular 
weight  distribution  of  solvent  chains,  1^1^  = 
11 8.1  kDa  ,  w®c  =0.7  and  IW®C2=287kDa,  w|c  = 
03. 

The  solid  curves  represent  the  continuous  GEX  for  the 
nonreactive  PDMS  soluble  chains  (m-PDMS)  whose 
parameters  are  listed  in  Table  1.  (B)  Bidisperse  approxi¬ 
mation  of  the  molecular  weight  distribution  of  network 
chains,  M^cl  =84.6kDa ,  w"c  =0.87  and  M^2  = 
260.2kDa .  i/VjC  =0.13.  The  solid  curves  represent  the 
continuous  GEX  for  the  reactive  network  precursor 
chains  (v-PDMS)  whose  parameters  are  listed  in  Table  1. 
[Color  figure  can  be  viewed  in  the  online  issue,  which  is 
available  at  wileyonlinelibrary.com.] 
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agreement  with  the  relaxation  processes  observed  in  the  exper¬ 
imental  dynamic  modulus  data  for  the  swollen  networks.  We 
find  that  dangling  chains  with  a  molecular  weight  of  187.5 
kDa  represent  the  experimental  G  data  reasonably  well. 


Fraction  of  dangling  chains 

Additionally,  the  volume  fraction  of  precursor  chains  that 
are  not  fully  cross-linked  during  the  curing  process  is  also 
unknown.  To  estimate  this  fraction,  we  use  results  of  Jensen 
and  coworkers.20  In  that  work,  it  was  shown  that  the  low- 
frequency  plateau  is  a  linear  function  of  the  density  of  cross¬ 
links  plus  trapped  entanglements.  Hence,  its  value  is  a  strong 
function  of  the  mass  fraction  of  active  structures  in  the  net¬ 
work,  in  agreement  with  classical  theories.  Using  the  DSM, 
several  systems  with  different  cross-link  densities  were  simu¬ 
lated,  and  an  empirical  relation  between  the  elastic  plateau 
and  the  fraction  of  the  dangling  structures  was  developed. 
In  this  work,  we  generalize  this  equation  to  a  bidisperse  net¬ 
work  as 


fro, ana 


WNC 
1—  W’sc 
pRT 


Hpr  <°-7«r.'- 


■IK 


,NC 


+ 


M"C2 


(0.74Zt,2  +  1K 


,NC 


(9) 


where  ZTi)l  is  the  average  number  of  trapped  entanglements, 
defined  as 


Ztj  — 


wnc 
1  —  Wsc 


(10) 


In  Eqs.  9  and  10,  w^c  is  the  fraction  of  network  chains, 
and  wsc  is  the  fraction  of  solvent  chains  which  is  known  as 
it  is  specified  for  a  given  loading  formulation.  For  the  calcu¬ 
lations  performed  in  this  article,  wsc  =0.5  and  vi’^c  and  w^c 
are  the  fractions  of  network  strands  of  the  bidisperse  distri¬ 
bution  of  network  strands  shown  in  Figure  2B. 

To  estimate  w’nc  ,  we  use  the  unloaded  network  elastic 
modulus  (Go=G'(cu  — >  0)  «  68kPa)  reported  by  Mrozek 
et  al.2  Using  Eq.  9,  our  estimate  for  /)  (whose  estimation  we 
discuss  in  Determination  of  the  DSM  Parameters  for  PDMS 
Section),  and  the  molecular  weight  distribution  discussed  in 
Blend  Composition  Section,  we  solve  Eq.  9  for  m’nc  •  For  the 
specific  system  being  considered,  the  fraction  of  network 
chains  is  found  to  be  wnc  =0.32  and  the  fraction  of  dan¬ 
gling  network  chains  is,  therefore,  w’dc  —  0.18. 


Determination  of  the  DSM  Parameters  for  PDMS 

The  three  input  parameters  for  the  DSM  are  the  number 
of  Kuhn  steps,  (Vk;  the  entanglement  activity,  /);  and  the 
characteristic  Kuhn  step  time,  tk-  Nk  can  be  calculated  using 
Nk=Mw/Mk,  where  MK  is  the  molecular  weight  of  a  Kuhn 
step,  reported  to  be  381  Da  for  PDMS.25 

The  DSM  parameters  for  PDMS,  ft  and  tk,  are  determined 
by  fitting  the  simulation  results  to  dynamic  modulus  experi¬ 
mental  measurements  for  the  pure  solvent  used  to  load  the 
network.  Dimensionless  relaxation  modulus  results  using  the 
DSM  are  shown  in  Figure  3A  for  each  of  the  two  molecular 
weights  in  the  bidisperse  molecular  weight  distribution  of 
the  solvent  chains.  A  modified  Baumgaertel,  Schausberger, 
and  Winter  (BSW)  function26  is  used  to  fit  the  simulation 
results  in  the  time  domain.  This  function  can  then  be  analyti¬ 
cally  transformed  to  the  frequency  domain  to  obtain  the 


Figure  3.  (A)  Relaxation  moduli  for  the  two  compo¬ 
nents  of  the  bidisperse  approximation  of  the 
entangled  solvent  (A/c  SCi  =16  and  WcS C2=43). 

The  symbols  represent  predictions  by  the  DSM  of  the 
individual  components  of  the  bidisperse  approximation 
of  the  molecular  weight  distribution  of  the  m-PDMS 
chains.  The  solid  lines  represent  a  fit  of  the  simulation 
results  by  a  BSW  spectrum.  (B)  Comparison  of  the  DSM 
prediction  for  the  pure  entangled  solvent  with  the  experi¬ 
mental  data  (symbols).  [Color  figure  can  be  viewed  in 
the  online  issue,  which  is  available  at  wileyonlinelibrary. 
com.] 

dynamic  modulus  shown  in  Figure  3B.  ft  was  estimated  to 
be  34.8.  As  the  entanglement  plateau  is  not  observed  in  the 
frequency  data,  there  is  some  uncertainty  in  this  value  but 
the  prediction  agrees  well  in  the  region  of  frequencies 
reported.  tk  was  simultaneously  estimated  during  this  fit  to 
experimental  data  to  be  0.026  ps.  Both  of  these  values  are 
architecture  and  molecular-weight  independent. 

Furthermore,  in  a  parallel  work  by  Andreev  et  al.,27  it  has 
been  shown  that  a  less  detailed  level  of  description  than  the 
DSM  can  be  used  to  get  quantitative  predictions  at  greatly 
reduced  cost.  The  Clustered  Fixed  Slip-Link  Model  is  the 
DSM  with  fixed  j8=l  and  Nc  instead  of  (VK  and  tc  instead  of 
tk-  Nc  as  a  function  of  N&  and  ft  is  given  by 

Mc  «  0.56(j6+1)Mk  (11) 

and  tc  as  a  function  of  tk  is  given  by 

tc  «  0.265y88/3TK  (12) 

Predictions  with  jS=  1  were  done  with  coarse-grained 
parameters  Nc>  nci  =  12  and  iVc,NC2  =  34,Ac,sci  =  16  and 
Mc,sc2=37,  and  (Vc,dc  =  25.  The  results  obtained  are  identical 
to  the  original,  unsealed  model. 
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Dynamic  Modulus  Predictions  of  the  Swollen 
Network 

Using  the  model  parameters  found  in  Determination  of  the 
DSM  Parameters  for  PDMS  Section,  we  can  now  construct 
the  self-consistent  CD  mean  field  according  to  Eq.  4.  Then, 
each  type  of  chain  in  the  ensemble  (irrespective  of  architec¬ 
ture  and  molecular  weight)  can  be  simulated  independently 
using  this  blend  pCD  field,  Eq.  4. 

Blend  predictions 

Figure  4A  shows  the  dimensionless  relaxation  moduli 
obtained  from  these  same  simulations  for  the  solvent  chains 
in  the  presence  of  the  blend  CD  field.  It  can  be  observed 
that  the  longer  molecular  weight  in  the  bidisperse  approxi¬ 
mation  has  a  long  relaxation  time  as  expected.  Figure  4B 
shows  the  dimensionless  relaxation  moduli  for  the  network 
chains  in  the  same  blend  CD  field.  The  relaxation  dynamics 
are  the  same  for  both  molecular  weights  as  in  the  network 
chains,  relaxation  occurs  only  from  CD,  and  these  CD  are 
the  same  for  both  molecular  weight  chains.  A  BSW  function 
is  used  to  fit  the  simulation  results  in  the  time  domain.  These 
functions  are  analytically  transformed  to  the  frequency 
domain  to  obtain  the  dynamic  modulus. 

Figure  4A  also  shows  the  dimensionless  relaxation  modu¬ 
lus  for  the  dangling  chains  in  the  presence  of  the  same  blend 
CD  field.  Compared  to  the  solvent,  the  dangling  chains  have 
a  much  longer  relaxation  time  as  they  are  fixed  at  one  end  to 
the  network,  and  SD  can  only  create  and  destroy  entangle¬ 
ments  at  one  end  of  the  chain. 

The  relaxation  modulus  for  each  architecture,  solvent 
chains,  dangling  network  chains,  and  network  chains,  can 
then  be  constructed  by  performing  a  weighted  average,  Eq. 
7,  of  the  relaxation  moduli  of  all  the  molecular  weights,  for 
each  architecture,  using  the  weight  fractions  found  in  Blend 
Composition  Section.  The  dynamic  modulus  is  then  obtained 
by  transforming  the  relaxation  modulus  to  the  frequency 
domain  analytically.  The  resulting  dynamic  moduli  for  each 
of  the  different  architectures  in  the  swollen  network  blend 
are  shown  in  Figure  5A.  Note  that  the  crossover  for  the  dan¬ 
gling  chains,  which  is  a  rough  indicator  of  the  longest  relax¬ 
ation  time,  occurs  at  much  lower  frequencies  than  that  for 
the  solvent  chains.  This  is  to  be  expected  as  the  dangling 
chains  are  similar  to  star-branched  melts,  which  can  not 
relax  by  SD  at  the  end  fixed  to  the  network.  For  the  dangling 
and  solvent  chains,  the  entanglement  plateau  can  be  observed 
at  high  frequencies.  It  can  also  be  observed  that  for  the  net¬ 
work  component  the  storage  modulus  and  the  loss  modulus 
do  not  crossover,  as  expected.  For  the  network,  an  entangle¬ 
ment  plateau  is  observed  at  high  frequencies,  and  an  elastic 
plateau  is  observed  at  low  frequencies.  Note  that  the  contri¬ 
bution  of  the  entangled  solvent  to  the  dynamic  modulus  of 
the  blend  is  only  relevant  at  high  frequencies  since  at  lower 
frequencies  the  solvent  chains  have  relaxed  significantly  and 
the  relative  magnitude  of  their  modulus  with  respect  to  net¬ 
work  and  network  dangling  chains  becomes  very  small. 
Observe  that  at  frequencies  around  101  Hz  the  G"  of  the 
swollen  network  blend  becomes  dominated  by  the  G"  of  the 
entangled  solvent.  The  overall  effect  of  the  high  molecular 
weight  solvent  is  to  increase  energy  dissipation  at  high  fre¬ 
quencies  while  having  a  negligible  effect  on  the  mechanical 
properties  at  lower  frequencies. 

The  dynamic  modulus  of  the  swollen  network  blend  is  cal¬ 
culated  using  Eq.  7,  with  h’sc  =0.5,  w^c  =0.18,  and 


Figure  4.  (A)  Relaxation  moduli  for  the  two  compo¬ 
nents  of  the  bidisperse  approximation  of  the 
SC,  and  of  the  monodisperse  approximation 
of  the  dangling  chains  (DC),  with  blend  CD 
Eq.  4. 

(B)  Relaxation  moduli  for  the  two  components  of  the 
bidisperse  approximation  of  the  network  chains  with 
self-consistent  blend  CD  Eq.  4.  [Color  figure  can  be 
viewed  in  the  online  issue,  which  is  available  at  wileyon- 
linelibrary.com.] 

w’nc  =0.32,  and  the  dynamic  moduli  for  each  of  the  individual 
architectures  shown  in  Figure  5A.  The  result  is  shown  in  Fig¬ 
ure  5B.  Recall  that  the  modulus  for  each  architecture  is  simu¬ 
lated  independently,  but  with  a  self-consistent  pCD  field  that 
accounts  for  the  presence  of  all  the  components  of  the  blend. 
The  agreement  with  the  experimental  data  is  good  and  could 
potentially  be  improved  by  increasing  the  number  of  molecu¬ 
lar  weights  in  the  discrete  approximation  of  the  GEX  distribu¬ 
tions  of  the  molecular  weights  of  each  of  the  components  of 
the  blend,  at  increased  computational  cost.  Noting  the  offset 
from  the  experimental  data,  suggests  that  our  estimate  for  fi 
may  be  a  bit  high.  Nonetheless,  we  observe  that  the  shape  is 
captured  well.  For  instance,  in  G' ,  the  slight  upturn  at  higher 
frequencies,  which  is  caused  by  the  increasing  contribution  of 
entangled  solvent  chains,  is  captured  correctly  by  the  model. 
In  G",  the  characteristic  power-law  behavior  at  intermediate 
and  low  frequencies  is  also  predicted  correctly  by  the  model. 

Predictions  of  Simple  Elongation  of  Swollen 
Network 

We  can  now  predict  the  nonlinear  rheology  without 
parameter  adjustments.  Simple  elongation  experiments  for 
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Figure  5.  Dynamic  modulus  of  the  swollen  network 
blend. 

(A)  Dynamic  modulus  of  each  component  in  the  blend 
with  self  consistent  blend  CD.  SC:  Solvent  chains;  DC; 
NC.  (B)  Comparison  of  the  DSM  prediction  for  the 
swollen  network  blend  with  the  experimental  data  of 
Mrozek  et  al.2  [Color  figure  can  be  viewed  in  the  online 
issue,  which  is  available  at  wileyonlinelibrary.com.] 


the  same  swollen  PDMS  network  were  reported  by  Sliozberg 
et  al.3  In  these  experiments,  a  dog-bone  shaped  sample  of 
initial  length  L  =  1  in.  is  elongated  at  cross-head  speeds 
between  0.05  and  5  in./s  using  a  MTS  810  servo-hydraulic 
load  frame.  In  this  type  of  experiment,  the  elongation  ratio 
of  the  sample,  X(X=l/L),  is  given  by 

X=Xt+ 1  (13) 

where  /  is  the  length  of  the  sample  during  the  experiment 
and  1  is  the  elongation  rate. 

The  elongation  ratio  is  related  to  the  Hencky  strain,  e, 

by22 

e=log  (A)  (14) 

and  consequently  the  Hencky  strain  rate  can  be  expressed  as 

A  i 

e=T - 7  =  (15) 

Xt+ 1  f+ 1/2 

which  is  used  in  the  DSM  algorithm  to  construct  the  velocity 
gradient,22  for  simple  elongation.21 

Here,  we  attempt  to  predict  the  experimental  results  for 
three  stretch  rates,  5,  0.5,  and  0.05  1/s.  The  simulation 
results  for  each  of  the  architectures  of  the  swollen  network 
blend  in  simple  elongation  for  the  highest  rate  are  shown  in 


Figure  6A.  Just  as  was  done  for  the  Green-Kubo  simulations 
presented  in  Dynamic  Modulus  Predictions  of  the  Swollen 
Network  Section,  in  this  case  too,  each  chain  in  the  ensem¬ 
ble  is  simulated  independently  in  a  self-consistent  CD  mean 
field,  Eq.  4,  that  accounts  for  the  presence  of  chains  with  dif¬ 
ferent  architectures  and  molecular  weights.  It  can  be 
observed  that  the  stress  contribution  from  network  chains  is 
much  larger  than  the  contribution  from  dangling  network 
chains  and  solvent  chains. 

The  average  swollen  network  response  is  obtained  by  cal¬ 
culating  the  weighted  average  of  the  individual  stresses 
shown  in  Figure  6A,  using  the  weight  fractions  found  in 
Blend  Composition  Section,  wsc  =0.5,  vt’oc  =0.18,  and 
w’nc  =0.32.  The  result  is  shown  in  Figure  6B,  together  with 
the  experimental  data  reported  by  Sliozberg  et  al.3  It  can  be 
observed  that  the  model  predicts  the  experimental  data  well 
with  small  deviations  for  1&5.  To  explain  the  fact  that  the 
high-molecular-weight  solvent  does  not  have  a  significant 
contribution  to  the  overall  mechanical  response  of  the  swol¬ 
len  gel,  we  estimated  the  Weissenberg  number  for  this  com¬ 
ponent  during  the  experiment  done  at  5  1/s  elongation  rate. 
For  the  solvent  component  of  the  system,  the  longest 


Figure  6.  (A)  DSM  prediction  for  simple  elongation  for 
each  of  the  components  of  the  PDMS  swol¬ 
len  network  blend  in  self-consistent  blend 
CD. 

The  figure  shows  the  elongational  stress,  a ,  as  a  function 
of  the  elongation  ratio,  L  (B)  Comparison  with  experi¬ 
mental  data  where  the  inset  shows  the  Weissenberg 
number  for  uniaxial  extension  for  the  solvent  component 
of  the  system  modeled  in  Determination  of  the  DSM 
Parameters  for  PDMS  Section,  1=5  1/s,  td=0.02s. 
[Color  figure  can  be  viewed  in  the  online  issue,  which  is 
available  at  wileyonlinelibrary.com.] 
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Figure  7.  Comparison  of  the  DSM  prediction  with 
experimental  data  for  elongation  rates 
1=5,  0.5,  and  0.05  1/s. 

The  symbols  are  experimental  data  and  the  solid  lines 
are  the  DSM  predictions.  [Color  figure  can  be  viewed  in 
the  online  issue,  which  is  available  at  wileyonlinelibrary. 
com.] 


relaxation  time  can  be  estimated  to  be  t,/  ~  0.02s  from  the 
cross-over  between  G'  and  G"  in  Figure  3B.  In  the  presence 
of  network,  this  longest  relaxation  time  goes  up  roughly  by  a 
factor  of  two  due  to  the  suppressed  CD  from  the  network 
chains,  as  can  be  observed  in  Figure  5A.  The  resulting  Weis- 
senberg  number  during  the  elongation  experiment,  Wi=e ij, 
is  shown  in  the  inset  of  Figure  6B.  Note  that  according  to 
this  analysis,  an  elongation  rate  on  the  order  of  100  1/s  or 
greater  would  be  required  to  observe  a  significant  effect  of 
the  high-molecular-weight  solvent  in  the  nonlinear  mechani¬ 
cal  response  of  the  swollen  network  blend.  Even  though  the 
direct  contribution  of  the  dangling  and  solvent  chains  is  very 
small,  there  is  an  indirect  contribution  via  CD  on  the  net¬ 
work  chains.  The  black  line  in  Figure  6A  includes  these 
effects. 

The  DSM  predictions  for  the  other  two  lower  elongation 
rates,  1=0.5,  and  0.05  1/s  tire  compared  with  experimental 
data3  in  Figure  7.  For  the  lowest  elongation  rate,  the  agree¬ 
ment  between  experiments  and  predictions  is  very  good.  How¬ 
ever,  there  is  a  significant  discrepancy  between  the  DSM 
prediction  and  the  experimental  data  for  the  intermediate  elon¬ 
gation  rate.  The  data  for  this  intermediate  elongation  rate  (0.5 
1/s)  does  not  appear  to  be  consistent  with  the  other  two  rates. 
The  stress  for  the  intermediate  rate  is  expected  to  lie  between 
the  stress  for  the  lowest  and  the  highest  rate.  Additionally,  the 
behavior  for  this  intermediate  rate  looks  more  Hookean  (lin¬ 
ear)  than  for  the  lowest  rate  which  also  indicates  that  the  data¬ 
set  for  0.5  1/s  seems  inconsistent.  The  inconsistencies  in  the 
data  may  have  been  caused  by  fracture  of  the  test  specimens 
(dog  bones)  during  the  elongation  experiments,  inhomogene¬ 
ities,  and  so  forth.  Regardless  of  these  inconsistencies  at  the 
intermediate  rate,  the  good  agreement  for  the  highest  and  low¬ 
est  rates  allows  us  to  conclude  that  the  DSM  can  successfully 
predict  the  nonlinear  elongation  behavior  of  PDMS  gels  swol¬ 
len  with  high  molecular  weight  solvent. 

Conclusions 

The  composition  of  gels  made  of  cross-linked  PDMS 
swollen  with  high-molecular-weight  PDMS  solvent  can  be 


tuned  to  produce  materials  with  a  large  storage  modulus  and 
an  enhanced  loss  modulus  at  high  frequencies.  More  specifi¬ 
cally,  the  mechanical  properties  of  swollen  network  blends 
depend  on  the  concentration  of  active  network  chains,  the 
concentration  of  solvent,  and  the  molecular  weight  distribu¬ 
tion  of  each  of  the  components.  To  characterize  such  a  sys¬ 
tem,  one  needs  a  theory  sufficiently  robust  that  it  can  be 
applied  to  linear  chains,  star-branched  chains,  and  cross- 
linked  gels  in  both  linear  and  nonlinear  rheology.  Although 
this  is  straightforward  in  the  DSM,  there  is  no  such  possibil¬ 
ity  in  tube  models.6'7  Moreover,  the  parameters  (j6,  tk)  can 
be  determined  from  dynamic  modulus  data  for  the  pure  sol¬ 
vent.  The  fraction  of  active  network  strands  was  estimated 
by  generalizing  previous  work  by  Khaliullin  and  coworkers20 
to  account  for  polydispersity  of  the  network  chains.  These 
parameters  were  then  used  to  predict  the  linear  and  nonlinear 
rheology  of  the  swollen  network. 

The  DSM  prediction  for  the  dynamic  modulus  of  the 
PDMS  swollen  gel  was  compared  to  experimental  data 
reported  by  Mrozek  et  al.2  We  find  that  the  DSM  produces 
good  quantitative  predictions  of  the  dynamic  modulus  of  this 
complex  system.  Our  results  show  that  the  effect  of  high- 
molecular-weight  solvent  becomes  relevant  in  the  overall 
dynamic  modulus  of  the  swollen  network  at  frequencies 
above  1  Hz.  In  this  frequency  range,  the  amount  of  energy 
dissipation  (G")  is  significantly  increased  due  to  the  presence 
of  the  entangled  solvent.  We  test  the  nonlinear  rheological 
predictions  of  the  DSM  against  elongation  experiments  of 
the  swollen  PDMS  networks.  We  find  that  the  DSM  can  pre¬ 
dict  well  the  experimental  observations.  Moreover,  valuable 
insight  into  the  specific  dynamics  of  each  component  in  this 
complex  blend  is  gained.  This  allows  us  to  conclude  that 
even  at  the  highest  elongation  rates  at  which  experiments 
have  been  performed,  the  contribution  of  the  high  molecular 
weight  solvent  to  the  overall  mechanical  response  of  the 
blend  is  negligible.  More  experimental  results  are  required  to 
validate  this  conclusion. 

The  good  quantitative  agreement  of  the  DSM  predictions 
with  experimental  results  indicates  that  the  model  can  be 
used  as  a  computational  tool28  to  design  swollen  gel  formu¬ 
lations  with  tailored  mechanical  properties  for  specific  appli¬ 
cations.  The  same  type  of  modeling  presented  in  this  work 
can  be  easily  further  developed  to  include  other  polymer 
architectures  such  as  star-branched  solvents,  which  will 
broaden  the  spectrum  of  mechanical  properties  that  can  be 
achieved  with  these  types  of  swollen  polymeric  gels. 
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